Quantum nondemolition measurements in a Paul 

trap 

A. Camacho * 
Astrophysikalisches Institut Potsdam. 
An der Sternwarte 16, D-14482 Potsdam, Germany. 



Abstract 

In this work a family of quantum nondemolition variables for the case of a 
particle caught in a Paul trap is obtained. Afterwards, in the context of the so 
called restricted path integral formalism, a continuous measuring process for 
this family of parameters is considered, and then the corresponding propaga- 
tors are calculated. In other words, the time evolution of a particle in a Paul 
trap, when the corresponding quantum nondemolition parameter is being con- 
tinuously monitored, is deduced. The probabilities associated with the possible 
measurement outputs are also obtained, and in this way new theoretical results 
emerge, which could allow us to confront the predictions of this restricted path 
integral formalism with the readouts of some future experiments. 



1 Introduction 



One of the more long-standing conundrums in modern physics is the so called quan- 
tum measurement problem (QMP) [1], the one from the very outset of quantum 
theory (QT) has provoked a deep interest by its unusual and paradoxical characte- 
ristics. Nowadays the attempts in the quest for a solution of this conceptual diffi- 
culty have been spawned by theoretical and also by practical reasons. On practical 
grounds, for instance, the necessity of detecting very small displacements in the case 
of gravitational-wave antenna [2], or the attempts to obtain high sensitivity in para- 
metric transducers [3] (a topic closely related to the design of gravitational radiation 
detectors), requires the analysis of QMP. 
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On the theoretical side, QMP plays a fundamental role in the understanding of 
the foundations of QT [1] , and therefore a better comprehension of this issue is closely 
related to the advancement of QT. 

Concerning QMP it is noteworthy to comment that one of the most interesting 
topics in this field is the so called quantum nondemolition measuring process [4], 
in which certain class of observables may be measured repeatedly with arbitrary 
precision. Of course, the dynamical evolution of the corresponding system limits this 
class of observables. Clearly this restriction is a direct consequence of the unavoidable 
back reaction of the measuring device upon the measured system. The fundamental 
idea behind a quantum nondemolition measurement (QNDM) is to monitor a variable 
such that the unavoidable disturbance of the conjugate one does not perturb the time 
evolution of the chosen parameter [3]. 

On experimental grounds QNDM is very promising [5], for instance, in the dy- 
namics of the interaction between measuring device and a mechanical oscillator the 
most important hurdles, that currently impede the achievement of the corresponding 
quantum regime, have already been identified, and the exploration of the quantum 
behavior in the context of macroscopic mechanical oscillators could bring in the near 
future decisive results [6]. 

In order to confront the theoretical predictions of RPIF against experimental 
outputs the necessity of performing continuous measurements upon a quantum system 
is imperative. In this respect we must mention that the current technology allows 
us to carry out repetitive measurements on a single quantum system. For instance, 
it is already possible [7] to confine and observe an individual electron in a Penning 
trap, or we can now also have a trapped single atom and observe its interaction with 
a radiation field, for example, by means of a laser fluorescence [8]. In the case of the 
so called Paul trap, which has led to the construction of a mass spectrometer [9], a 
ion is trapped employing a high-frequency electric quadrupole field [10]. As is known, 
this idea can also be extended to the case of neutral atoms, laser cooled and stopped 
atoms are confined in a magnetic quadrupole trap formed by two opposed, separated, 
coaxial current loops [11]. 

The aforementioned experiments open the possibility of confronting, in a near 
future, the theoretical predictions of some formalisms that claim to describe the 
interaction between measuring apparatus and measured system against experimental 
results. One of these formalisms is the so called Restricted Path-Integral Formalism 
(RPIF) [12]. The main idea in this approach is the restriction by means of a weight 
functional of the integration domain of the path-integral that renders the corres- 
ponding propagator of the analyzed system, when one or more of its parameters are 
subject to a continuous measurement process. 

Let us explain this point a little bit better, and suppose that we have a particle 
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which has a one-dimensional movement. The amplitude U(q",q') for this particle 
to move from the point q' to the point q" is called propagator, and it is given by 
Feynman [13] 

U(q",q') = J d[q}exp Q-%]) , (1) 

here we must integrate over all the possible trajectories q(t), and S[q] is the action of 
the system, which is defined as 

S[q]= f dtL(q,q). (2) 
Jt' 

Let us now suppose that we perform a continuous measurement of the position 
of this particle, such that we obtain as result of this measurement process a certain 
output a(t). In other words, the measurement process gives the value a(t) for the 
coordinate q(t) at each time t, and this output has associated a certain error Aa, 
which is determined by the experimental resolution of the measuring device. The 
amplitude U[ a ](q",q') can be now thought of as a probability amplitude for the con- 
tinuous measurement process to give the result a(t). Taking the square modulus of 
this amplitude allows us to find the probability density for different measurement 
outputs. 

Clearly, the integration in the Feynman path-integral should be restricted to 
those trajectories that match with the experimental output. RPIF states that this 
condition can be introduced by means of a weight functional u)[ a ] [14]. This means 
that expression (1) becomes now under a continuous measurement process 

U[ a ] = J d[q]u>[ a ]exp(iS[q]). (3) 

The more probable the trajectory [q] is, according to the output a, the bigger 
that Cow,), becomes [12], this means that the value of a>[ a i is approximately one for all 
trajectories, [q], that agree with the measurement output a and it is almost for 
those that do not match with the result of the experiment. 

Clearly, the weight functional contains all the information about the interaction 
between measuring device and measured system, and the problems that in this context 
appear are two: 

(i) The concrete form of the weight functional oj[ a ] depends on the measuring 
device [15], in other words, the involved experimental constructions determine these 
weight functionals. We must determine uj[ a j starting from the knowledge that we have 
of the measuring device, a non-trivial problem. 
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(ii) If we wish an analytical expression for L/[ a ] , then the resulting weight functional 
must render an analytically handleable functional integral. 

In this work we will consider the case of a particle caught in a Paul trap, and 
the nonlinear differential equation, the one defines the corresponding quantum non- 
demolition variables, will be solved. As solution a family of quantum nondemolition 
parameters will be obtained. Afterwards, a continuous measuring process for one of 
these parameters will be considered, and the associated propagators will be calcu- 
lated. Finally, the probabilities of the different measurement outputs are obtained. 



2 Quantum nondemolition variables in a Paul trap 



It is already known that in the case of a particle in a Paul trap we have a harmonic 
oscillator, the one possesses a frequency equal to U — Vcos(ujt), being U, V, and uj 
constants which depend on the electric quadrupole field used to trap the particle. 
The motion equations of an electrically charged particle caught in a Paul trap are 
[10] (here we will assume that ions are injected in the y-direction and that there is 
electric field only along the x- and ^-coordinates) 

+ - Vcos(iot)]x(t) = 0, (4) 

z(t) - - Vcos(ujt)}z(t) = 0. (5) 

here e is the electric charge of the particle, m its mass, and 2r the distance between 
the electrodes that conform part of the experimental apparatus. The solutions to 
these equations are the so called Mathieu functions [16]. 

Let us now first consider the motion only along the X clXIS. This is no restriction 
at all, because the complete motion can be separated in two independent motions. 
Thus, our starting point is the following Lagrangian 

1 1 

L = -mx 2 (t) - -m[U - Vcos(ujt)]x 2 (t). (6) 

here we have introduced the following definitions: U = -^U and V = -%V\ 

In the case of a harmonic oscillator with a time dependent frequency a quantum 
nondemolition (QND) variable may be found in the form [12] 
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A = pp + aq, 



(7) 



where p and q denote the momentum and position, respectively, and p and a satisfy 
the following differential equation [12] 



d L = !l 

dt m 



+ m[U - Vcos(ujt)\, 



(8) 



here / 



Let us now consider the following function: 

„ in dx 



x(t) dt 

where x(t) is any of the possible solutions to 



(9) 



x(t) + [U- Vcos(ut)]x(t) = 0. (10) 

It is straightforward to check that F is a solution to (8). Hence, with the solutions 
of the motion equation we may easily find QND-variables. In other words, we have 
found a family of QND-parameters for the case of a particle in a Paul trap. 



A(t) = a(t) 



1 dx 



— m- 



-q + p 



(11) 



x(t) dt 

If we set V = 0, then we recover, from expression (9), the situation already 
known of a harmonic oscillator with time independent frequency [12]. Indeed, in that 
case, x(t) = —cos(VUt), and then p/a = my/17 tan(v / f/t), where VH = u is the 
corresponding frequency of oscillation. 



Continuous measurements, propagators, and pro- 
bability densities 



Let us now consider the situation in which our QND-variable, A, is being conti- 
nuously monitored, here we set a(t) = 1. According to RPIF [12] the corresponding 
propagator (when the measurement output reads a(t)) is given by 
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U [a] = jP d[q}d[p] exp ji £ Q?l + _L[tf - Vcosiujt)^ <ft J u) [a] . (12) 

As was mentioned at the end of the first section, the whole interaction between 
measuring apparatus and measured system is contained in Cj^ . At this point we must 
choose a particular weight functional, and it will be a gaussian one 

£[«] = ^p {" /' [^(*) " «(f)] 2 *} • (13) 

The reasons for this choice lie on the fact that the results coming from a Heaveside 
weight functional [14] and those coming from a gaussian one [17] coincide up to the 
order of magnitude. These last remarks allow us to consider a gaussian weight func- 
tional as an approximation of the correct expression. Hence, it will be supposed that 
the weight functional of our measuring device has precisely this gaussian form. We 
may wonder if this is not an unphysical assumption, and in favor of this argument we 
may comment that recently it has been proved that there are measuring apparatuses 
which show this kind of behavior [18]. 

Therefore, expression (12) becomes 



fit, - /; 4144 - P { i /; + "-^f + ^[A - «?) *} .04, 

This last path integral is gaussian in p and q, and therefore, it can be easily 
calculated [19] 



^ = H~r^/* a2dt 



f TAa 2 - 2imk ft" (a^) 2 [2m 2 ha + imTAa 2 f3] , 

where we have introduced three new definitions, namely, a = (-) 2 + U — Vcos(u;t), 
P = U- Vcos(cot), and 7 = Am 2 h 2 + T 2 Aa 4 . 

The probability densities, associated with the different measurement outputs, read 
(according to the expression P[ a ] = |£7[ a ]| 2 ) 



6 



P [a] = exp 



TAa 2 Jt 



a dt 



exp 



TAa 2 (af) 2 (f ) 2 + 2/? 
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Jt' 



2 T 2 Aa 4 Q2 



a 2 + 



(16) 



4 Conclusions 



In this work we have considered the case of a particle caught in a Paul trap, and, 
after solving the corresponding nonlinear differential equation, a family of quantum 
nondemolition variables has been obtained, expression (11). Afterwards, a continuous 
quantum measurement for an element of this family was considered, and, along the 
ideas of the so called restricted path integral, the corresponding propagator has been 
calculated, expression (15). Finally, the associated probability densities were derived, 
expression (16). The present work complements a previous paper in which a quantum 
demolition measurement for a particle in a Paul trap was analyzed [20]. 

Clearly, we may notice in our last equation that there is no standard quantum 
limit, in other words, we may measure A with an arbitrarily small error, and in 
consequence all the necessary information can be extracted. 

Looking at (16) we may notice that in the limit Aa — > we obtain P[ a ] — > 0. 
In other words, in the limit of very precise measurements all the possible readouts 
have the same probability density. This last fact is a quantum feature, indeed, in the 
nonquantum case only the solution to the classical motion equations has a nonvan- 
ishing probability. This last remark does not mean that a strong disturbance of the 
corresponding observable is present. Indeed, we may find this kind of behavior even 
in a much simpler situation, namely, in the case of a harmonic oscillator the limit of 
very precise measurements (Aa — > 0) renders also the result Pr a i — > (see equation 
(6.32) in [12]). This characteristic is also present in the situation when a quantum 
nondemolition variable, for the case of a particle moving in an inhomogeneous gra- 
vitational field, is monitored in the limit of very small intrumental error [21]. The 
opposite case, namely, the limit of rough measurements, Aa — > oo, renders also the 
result P[ a ] — > 0. 

Of course, we could have a very small experimental error, but the case Aa — > 
is an idealization, and this limit has to be understood in the sense that if the expe- 
rimental resolution is much smaller than all the relevant physical variables, then we 
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could expect to have a probability independent of the measurement outputs. Expe- 
rimentally this case could be a very difficult one, consider, for example, the current 
experimental resolution in the case of Paul or Penning traps [10], which lies very far 
from this idealization. 

There are already some theoretical results [22, 23] that provide a framework which 
could allow us to confront the predictions of RPIF against some possible future expe- 
riments. Nevertheless, the present work renders new theoretical predictions coming 
from RPIF, which in the future could be confronted with experimental readouts. 
For example, expression (16) predicts a very particular dependence (on the involved 
measurement readouts) for the ratio of the probability densities associated with two 
different experimental results. Indeed, if b(t) ^ a(t), then its is readily seen that 
P[a]/P[b] = di - 2 — b 2 ), where g is a real function 



TAa' 



ft" / X 



TAa 2 ( -«"(a 2 -6 2 )(f)H(f) 2 + 2/3 



^ 7 I, a 2 + 



(17) 



To finish, let us comment an additional characteristic of expression (16). We know 
that there are some QNDMs in which an absolute limit may appear. For instance, if 
the linear momentum of a free particle is monitored, this aforementioned limit may 
emerge, when the instantaneous measurement of position of the test particle, before 
and after the monitoring of the linear momentum, is done (see page 99 of reference 
[12]). Clearly, position and linear momentum are canonical conjugate variables to 
each other, and that is why this absolute limit appears. At this point we may wonder 
why (16) shows no absolute limit, in our case, position, before and after monitoring 
of A(t), has also been instantaneously measured, i.e., q' and q" are present from the 
very begining in our mathematical expressions, see (12). 

The answer to this question stems from the fact that in our case the measured 
quantity, A = a(t)[—m-^^q+p], is not the canonical conjugate variable of position, 
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